A general theoretical description of N-body recombination 
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A formula for the cross section and event rate constant describing recombination of TV particles 
are derived in terms of general S-matrix elements. Our result immediately yields the generalized 
Wigner threshold scaling for the recombination of TV bosons. A semi-analytic formula encapsulates 
the overall scaling with energy and scattering length, as well as resonant modifications by the 
presence of TV-body states near the threshold collision energy in the entrance channel. We then 
apply our model to the case of four-boson recombination into an Efimov trimer and a free atom. 



Few-body processes have played an increasingly cru- 
cial role in the physics of strongly interacting quantum 
gases. Three-body recombination in particular [JJ con- 
tributes strongly to atom-loss and is primarily respon- 
sible for controlling the lifetime of Bose-Einstein con- 
densates (BEC) because the kinetic energy released in 
the reaction is usually sufficient to eject the collision 
partners from the trapping potential [2 . At ultracold 
temperatures, the physics of three-body recombination 
is largely controlled by the effective long-range potential 
in the entrance (three-body) channel. Hence, a threshold 
analysis of these potentials immediately yields informa- 
tion about the scaling behavior of the event rate constant 
K$ with respect to both the energy E and the (s-wave) 
scattering length a. For instance, a threshold analysis 
shows that K3 scales as a 4 and E° for bosons [TJ, mak- 
ing three-body recombination an important process in 
the ultracold limit. In contrast, the long-lived nature of 
the ultracold polarized Fermi gas is explained by the E 2 
scaling of K3 near threshold [3 . Further, the presence of 
Efimov states [4] near the threshold collision energy reso- 
nantly enhances recombination [1] . This resonant feature 
was exploited to confirm the existence of Efimov states 
in experiments with ultracold cesium [5]. The increas- 
ing number of experiments observing Efimov physics [6 
highlights the importance of few-body physics in our un- 
derstanding of strongly interacting quantum gases. 

The natural extension of three-body recombination is 
to four-bodies. Considerable progress has in fact been 
made towards the calculation of four-body processes, no- 
tably collision cross sections involving two-body fragmen- 
tation channels [7] and bound state energies [Sj. Scatter- 
ing processes involving four or more free atoms are far 
more complex and require a deep analysis of the mul- 
tiparticle continua [5J [TO]. Nevertheless, recent studies 
have proposed four-body recombination as an efficient 
process for the production of Efimov trimers [5J ITU] , 
Four-boson borromean states (four-body bound states 
with no bound subsystems) with universal properties tied 
to Efimov physics were recently predicted further 
explored |10j , and measured by the Innsbruck group |12j 



in atom loss due to recombination. The present analysis 
is concerned with the recombination of TV particles into 
a channel with TV — 1 particles bound together and one 
free. 

We use hyperspherical coordinates (see [2] and [TS] 
for useful reviews), in which Jacobi vectors pi are trans- 
formed to a set of angular coordinates collectively de- 
noted f2, plus a radial coordinate called the hyperra- 
dius R defined by ^nR 2 — YliLi 1 fHPii where /zjv = 
[(n i m t)/ M } 1/iN ~ 1) is the TV-body reduced mass, M = 
^2 , rrii is the total mass of the system and /ij is the re- 
duced mass associated with the i th Jacobi vector. At 
large R, the solutions to the angular portion of the 
Schrodinger equation yield the fragmentation channels of 
the TV-body system, and the quantum numbers labeling 
those solutions index the ^-matrix. 

Derivation of the generalized cross section. — This for- 
mulation begins by considering scattering by a purely 
hyperradial potential in d-dimensions, and then obtains 
the generalized cross section "by inspection" . For clarity, 
we adopt a notation that resembles the usual derivation 
in three dimensions. 

In d-dimensions, the wavefunction at large R behaves 
as: 

^ s pikR 
^^^ R + f(^') W=m (1) 

Equivalently, an expansion in hyperspherical harmonics 
is written in terms of unknown coefficients Ax^: 

m 11 = Y J A^Y^{R)[f x (kR) cosri A - n d x (kR) smS x ] 

X,/.t 

(2) 

Here, Y\ M are hyperspherical harmonics (solutions to the 
free-space angular equation (A 2 — A(A + d — 2))Y\^ = 0, 
where A 2 is the grand angular momentum operator |15j ) 
and jf (nf) are hyperspherical Bessel (Neumann) func- 
tions [TS]. 

Identification of the incoming wave parts of 'J 7 and "i 11 
yields the coefficients A\^, whose insertion into Eq. (p| 
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gives: 



S'-matrix elements as 



f(k,k') = 



2tt 



jA e -i(d/2-l+A)7r/2-47r/4 

A// 

xY^(k)Y Xfl (k')(e^-l). (3) 



The immediate generalization of this elastic scattering 
amplitude to an anisotropic short-range potential is: 



f(k,k') 



E i^-^Y^Yy^k') 

XfiX' fi' 
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where ?/a = (d/2 — 1 + \)ir/2 + ir/4. Upon integrating 
|/(fc, fc')| 2 over all final hyperangles fc, and averaging over 
all initial hyperangles k' as would be appropriate to a 
gas phase experiment, we obtain the average integrated 
elastic scattering cross section by a short-range potential: 
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where n(d) = 2ir d / 2 /T{d/2) is the total solid angle in d- 
dimensions [15]. This last expression is immediately in- 
terpreted as the generalized average cross section result- 
ing from a scattering event that takes an initial channel 
into a final channel, i = A'/i' — > Xfj, = /. Since this 
S-matrix is manifestly unitary in this representation, it 
immediately applies to inelastic collisions as well, includ- 
ing TV-body recombination. It is worth noting that the 
sum in Eq. (|5| should include degeneracies, and that the 
cross section should be appropriately averaged over ini- 
tial spin substates and collision energies. 

In this form, we can interpret the generalized cross 
section derived above in terms of the unitary S-matrix 
computed by solving the exact coupled-channels refor- 
mulation of the few-body problem within the adiabatic 
hyperspherical representation |16J . In principle this can 
describe collisions of an arbitrary number of particles. 
Identical particle symmetry is handled by summing over 
all indistinguishable amplitudes before taking the square, 
averaging over incident directions and momenta for a 
given energy, followed by integrating over distinguishable 



final states to obtain the total cross section: a 



indist 



N p a 



Here N p is the number of terms in the permu- 
tation symmetry projection operator (e.g. for TV identical 
particles, N p — TV!.) 

The cross section for total angular momentum J and 
parity II includes an explicit 2 J + 1 degeneracy. Hence, 
the total generalized cross section (with dimensions of 
length^ -1 ) for TV particles in all incoming channels i to 
scatter into the final state /, properly normalized for 
identical particle symmetry, is given in terms of general 
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The event rate constant [recombination probability per 
second for each distinguishable TV-group within a (unit 
volume)*^ -1 )] is the generalized cross section Eq. ^ 
multiplied by a factor of the iV-body hyperradial "veloc- 
ity" (including factors of Ti to explicitly show the units 
oiK N ): 
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Treatment of N bosons. — The 
a relevant S'-matrix element from an 
batic hyperspherical viewpoint is seen to 
Sj" = ($/ (R;Q.)F f (R)\S\$f(R; Sl)Fi(R)), where 

<&/ and $j are the channel functions (i.e. solutions 
to the hyperangular part of the Schrodinger equation 
in the limit R — > oo) for the entrance and final channel 
respectively. For all TV-body entrance channels we 
have i — * A, where A is the hyperangular momentum 
quantum number associated with cigcnfunction Y\ u . 
The functions Fi and Ff are the large R solutions 
to the coupled hyperradial equations obtained in the 
adiabatic hyperspherical representation [16 (in units 
where h = 1): 



-1 d 2 
2/^dR 2 



Wi(R)-E 



F i + J2Vif(R)F f =0, (8) 



Asymptotically (as R — > oo) the couplings Vif(R) van- 
ish and the effective potentials in the N-body continuum 
channels arc expressed in terms of an effective angular 
momentum quantum number l e : 



W X (R) 



h(le + 1) 



with l e = (2A + d-3)/2. (9) 



2/j, N R 2 

Near threshold, the recombination cross section is con- 
trolled by the lowest TV-body entrance channel i = A — > 
Amin = 0. For TV identical bosons in a thermal (non- 
quantum-degenerate) gas cloud, the dominant contribu- 
tion to Eq. ^ is from J u = + . Further, unitarity of 
the S-matrix allows us to write the total rate constant 
summed over all final channels (using d — 3TV — 3) as: 
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Based on prior hyperspherical treatments for three |19j 
and four [10] bosons, in cases where one or more (TV — 1)- 
body state is bound with no bound states of fewer par- 
ticles, we expect the hyperradial potentials to schemat- 
ically look like those sketched in Fig. [l] Monte Carlo 
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FIG. 1: (color online) A schematic representation of the 
iV-boson hyperradial potential curves is shown. When a 
metastable iV-boson state crosses the collision energy thresh- 
old at E = 0, iV-body recombination into a lower channel 
with N — 1 atoms bound plus one free atom is resonantly 
enhanced. 



calculations [T3] of A^-boson systems provide evidence of 
super-bori'omean states which could produce the topol- 
ogy shown in Fig. [T] For four bosons, potentials with 
the topology shown occur at negative scattering lengths 
(a < 0). 

The semiclassical (WKBJ) treatment of Berry [20] can 
be generalized, giving a complex phase shift Si c for colli- 
sions of N identical bosons. The phase shift is acquired 
in terms of the phase 4> accumulated in the classically al- 
lowed region (i?3 < R < R2), and the tunneling integral 
7 in the classically forbidden regions: 



q(R)dR; and 7 = Im 



(3JV-5)/2fe 



q(R)dR, 



(11) 

where q(R) = y/2n N (E -W'{R)) (W is the effective 
potential including the Langer correction [H]), and R* 
is the hyperradius at which the inelastic coupling to the 
exit channel peaks (typically at hyperradii much greater 
than the range of the two-body interaction ro; i.e. ro <C 
R* ^ R3 for recombination into a universal (N— l)-body 
state.) 

The inelastic amplitude is incorporated into the phase 
by letting (f> — > (f> + ir\ [22J. Applying the connection 
formulas for each of the classical turning points Ri, R2 
and i?3 shown in Fig. [I] gives the phase shift: 



Si. = S, 



Here, <5,^ is the real semiclassical phase shift derived by 
ignoring the interior region, which does not contribute 
to the inelastic probability. The S-matrix element de- 
scribing scattering from one A^-boson state to another is 
then simply: Sqo = e 2lSl ^ , and the total probability to 
recombinc into all available final channels is: 

1 10 12 e ~ 27 sinh(2?y) 
1-ISool =— — — — %T,W (13) 
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Qe~ 27 cot(0 + 7>7- tt/2)^ (12) 



where the function A(ri, 7, <j)): 
4(77,7,0) 



e" 27 

1 H — tanh (77 + icp) 



(14) 



is equal to unity in the threshold limit unless both 77 — ► 
and <fi — > 7r/2. 

The threshold energy dependence of is found by 
breaking 7 into two pieces corresponding to the regions 
shown in Fig. [T| 

p{3N-5)/2k ra\a\ 

ln= / \q{R)\dR; Ji = lm q(R)dR. 

Ja\a\ JR* 

(15) 

For recombination into a universal (N — l)-body bound 
state, we expect R* ^> ro, in which case 7, (f> and n 
are also universal. A simple calculation shows e~ 21 " — 
(2fccv|a|/(3A — 5)) 3Ar ~ 5 . It is convenient to introduce an 
a-independent function C(a) — C^a^ 1 ^'^ e^ 211 . The 
constant Cn must be adjusted to give the correct over- 
all scale of A'Jy . The full A^-boson recombination rate 
constant in the ultracold limit is then: 
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C(a) sinh (277) 

cos 2 <p + sinh 2 77 
(16) 

Note that this is a quantitative result valid in the thresh- 
old regime. It is a constant (independent of k), and scales 
roughly as |a| 3Ar_5 (in agreement with [IB])- Hence, N- 
body processes could contribute to the total atom loss 
at higher particle density n through terms of the form 
—K N n N : t| = Y^n^2 ~I J N'n N , where the atom loss 
rate Ln is related to the event rate by Ln = Kn /(N— 1)! 
provided each recombination event ejects A^ atoms from 
the trapped gas. (Recall also that for a quantum- 
degenerate Bose gas, the above expressions for K~n must 



be divided by AH [T7]-) For N = 3, Eq. (Jl6| is consistent 
with the expression found by Braaten and Hammer |22j 
and the phase (f> acquires the universal log-periodicity 
characteristic of Efimov physics 0]. 

Recombination of four bosons with a < 0. — Now we 
specialize our results to the case of four identical bosons 
(d = 9, J = 0, Ap = 4!, and no degeneracy.) Using nu- 
merical hyperradial potential curves for four bosons inter- 
acting via a short-range model potential [lOj , we obtain 
K\ both by solving the coupled channels numerically 
and by using Eq. ( 16 1 specialized to N = 4. In Fig. [2] 



cos 2 <p + sinh 77 



we show K® on an absolute scale [cm 9 /s] (where we 
use parameters appropriate for 133 Cs: m = 244188 a.u. 
and 7*0 = 100 a.u.), and in model units of [?irg//j 4 ]. The 
horziontal axis is shown in units of ro and in units of the 

universal "three-body parameter" k = ^2713 \E^g\/h, 
(2) 

where E^g is the bound-state energy of the second Efi- 
mov trimer at unitarity. The overall magnitude of K® 
is governed by — a/ro, while the position of the peaks 



4 



0.15 0.23 0.3 0.45 

io- 23 



0.75 1.1 




10 26 




10 23 








10 20 








10 17 




10 14 




10 11 





50 70 



FIG. 2: (color online) The four-boson recombination rate 
constant K% + is shown over approximately one range of the 
discrete scaling factor e 71 ^ 30 [10]. The dots are obtained by 
numerically solving the coupled hyperradial equations, while 
the solid green curve is obtained from Eq. ( 16 1 with N = 4, 
r\ — 0.01 and Ca = 55. The dashed red curve is a calculation 
using Eq. (161, but ignoring the additional suppression due 
to 7 /. T he inset shows the universal phases <f> and 7/ used in 
Eq. (16 1. While 7/ is shown here for a — 10, results for K\ 



are independent of a (see text). 



is fixed with respect to — na. However, because k and 
ro are related by a nonuniversal factor, the relation- 
ship between the two horizontal axes in Fig. [2] is model- 
dependent. These results show good agreement between 
the two calculations demonstrating the validity of the 
WKBJ approximation The results show the overall \a\ 7 
scaling modified by resonant peaks at: 



0.67 and 



1.40 



(17) 



A cusp in K4 appears at fc|a| ~ 1.56 when a new atom- 
trimer channel appears (i.e. a new Efimov state becomes 
bound), whereas peaks appear when a four-boson state 
sits at the threshold of the entrance channel. 

As was recently noted [It)], two four-boson states are 
bound at slightly less negative values of a than the values 
of a at which an Efimov trimer becomes bound. In the 
potentials of Fig. [T] for N = 4, the entrance channel is 
capable of supporting four-boson bound states, but the 
second trimer-atom channel is not yet available. Because 
inelastic transitions occur at hyperradii of order the size 
of the lowest Efimov trimer, C, <j>, and 77 are universal 
functions of a. We assume 77 is independent of a, and use 
7/ = 0.01 and C4 = 55 in Fig. [2] For larger values of |a|, 
the universal structure of Fig7|2] repeats with the three- 
boson scale factor e 7 ^ ss 22.7 [ID]. The Innsbruck group 
has recently observed this universal connection between 
Efimov states and four-boson states and confirmed the 



spacing implied by Eq. (17 1, and the observed atom loss 
rates are consistent with the absolute scale in Fig. [2] 12J . 

In conclusion, we have derived a general formula for 
the event rate constant for ./V-body recombination. The 
generalized Wigner threshold scaling laws immediately 



follow, and the overall scaling is resonantly modified by 
the presence of metastable TV-body states near threshold. 
We then specialize to four-bosons with a < 0. 

This work was funded in part by the National Science 
Foundation. The authors thank B. D. Esry for extensive 
discussions related to the present derivations. 



* Electronic address: mehtan@iilaul.colorado.edu 
t Electronic address: jpdincao@jila.colorado.edu 

[1] B. D. Esry, C. H. Greene, and J. P. Burke Jr, Phys. Rev. 
Lett. 83, 1751 (1999); P. O. Fedichev, M. W. Reynolds, 
and G. V. Shlyapnikov, ibid. 77, 2921 (1996); E. Nielsen 
and J. H. Macek, ibid. 83, 1566 (1999); P. F. Bedaque, 
E. Braaten, and H.-W. Hammer, ibid. 85, 908 (2000); 
E. Braaten and H.-W. Hammer, ibid. 87, 160407 (2001). 

[2] E. Braaten and H. Hammer, Phys. Rep. 428, 259 (2006). 

[3] H. Suno, B. D. Esry, and C. H. Greene, Phys. Rev. Lett. 
90, 053202 (2003). 

[4] V. Efimov Phys. Lett. B 33, 563 (1970); 
Sov. J. Nuc. Phys. 10, 62 (1970). 

[5] T. Kraemer et al, et al., Nature 440, 315 (2006). 

[6] S. Knoop et al, Nature Physics 5, 227 (2009); T. B. 
Ottenstein et al, Phys. Rev. Lett. 101, 203202 (2008) 
J. H. Huckans et al, Phys. Rev. Lett. 102, 165302 (2009) 
G. Barontini et al, Phys. Rev. Lett. 103, 043201 (2009) 
M. Zaccanti et al, Nature Physics 5 586 (2009); 

[7] D. S. Petrov, C. Salomon, and G. V. Shlyapnikov, Phys. 
Rev. Lett. 93, 090404 (2004); N. P. Mehta et al, (2007), 
arXiv:0706.1296; J. P. D'Incao et al, Phys. Rev. A 79, 
030501 (R) (2009); A. Deltuva and A. C. Fonseca, Phys. 
Rev. C 75, 014005 (2007) J. P. D'Incao, J. von Stecher, 
and C. H. Greene Phys. Rev. Lett. 103, 033004 (2009). 

[8] L. Platter, H.-W. Hammer, and Ulf-G. Meifiner, Phys. 
Rev. A 70, 052101 (2004). 

[9] Y. Wang and B. D. Esry, Phys. Rev. Lett. 102, 133201 
(2009). 

[10] J. von Stecher, J. P. D'Incao, and C. H. Greene (Nature 

Physics 5, 417 2009) 
[11] H.-W. Hammer, L. Platter, Eur. Phys. J. A 32, 113 

(2007). 

[12] F. Ferlaino et al, Phys. Rev. Lett. 102, 140401 (2009). 
[13] G. J. Hanna and D. Blume, Phys. Rev. A 74, 063604 
(2006). 

[14] Y. F. Smirnov and K. V. Shitikova, Sov. J. Part. Nucl. 
8, 344 (1977). 

[15] J. Avery, Hyperspherical Harmonics: Applications in 
Quantum Theory (Kluwer Academic Publishers, Norwell 
MA, 1989), see pp. 52-53 for the d-dimensional plane- 
wave expansion and the definition of j*. 

[16] J. H. Macek, J. Phys. B 1, 831 (1968). 

[17] Y. Kagan, B. V. Svistunov, and G. V. Shlyapnikov, JETP 
Lett. 42, 209 (1985). 

[18] M. Th0gersen, D. V. Fedorov, and A. S. Jensen, Euro- 
phys. Lett. 83, 30012 (2008). 

[19] J. P. D'Incao and B. D. Esry, Phys. Rev. Lett 94, 213201 
(2005). 

[20] M. V. Berry, Proc. Phys. Soc. 88, 285 (1966). 

[21] P. M. Morse and H. A. Feshbach, Methods of Theoretical 

Physics (McGraw-Hill, New York, 1953), 1st ed. 
[22] E. Braaten and H.-W. Hammer, Phys. Rev. A 70, 042706 



(2004). 



